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We report the results of a search for the bottomonium ground state rib{lS) in the photon energy 
spectrum with a sample of (109 ± 1) million of T{3S) recorded at the T{3S) energy with the BABAR 
detector at the PEP-II B factory at SLAC. We observe a peak in the photon energy spectrum at 
-E7 — 921.2l2'8(s^^^) =t 2.4(syst) MeV with a significance of 10 standard deviations. We interpret 
the observed peak as being due to monochromatic photons from the radiative transition T{3S) —^ 
7»76(15'). This photon energy corresponds to an r;6(lS') mass of 9388.9I2 3(stat)±2.7(syst) MeV/c^. 
The hyperfine r{lS)-rit{lS) mass splitting is 71.4l3;i(stat) ± 2.7(syst) MeV/c^. The branching 
fraction for this radiative T(3S) decay is estimated to be (4.8 ± 0.5(stat) ± 1.2(syst)) x 10~*. 

PACS numbers: 13.20.Gd, 14.40.Gx, 14.65.Fy 



Thirty years after the discovery of the narrow T{nS) 
resonances T], no evidence has been reported for the 
spin-singlet pseudoscalar partners rib{nS) of these states. 
Measurement of the hyperfine mass splittings between 
the triplet and singlet states in quarkonium systems is 
of key importance in understanding the role of spin-spin 
interactions in quarkonium models and in testing QCD 
calculations [3] ■ Theoretical estimates of the mass split- 
ting between the IS* singlet and triplet states vary from 
36 MeV/c2 to 100 MeV/c^ |3|. 

In this letter, we report the observation of the radiative 
transition T{?>S) -^ jrjbilS), where the rib{lS), hereafter 
referred to as the rjb, is the pseudoscalar partner of the 
triplet state T{IS), and corresponds to the ground state 
of the bottomonium system. Theoretical predictions of 
the decay branching fraction range from 1 to 20 x 10"'* [3|, 
where the unknown rjb mass is a major source of the un- 
certainties. The current limit from the CLEO III exper- 
iment, B[r{3S) ^jr]b] < 4.3 x 10"^ at 90% confidence 
level, is based on 1.39 fb"^ of r(35) data [4]. 

The data sample used in this study was collected with 
the BaBAR detector [5| at the PEP-II asymmetric-energy 
e+e" storage rings. It consists of 28.0 fb""'^ of integrated 
luminosity collected at a e+e~ center-of-mass (CM) en- 
ergy of 10.355 GeV, corresponding to the mass of the 
T(3S) resonance. Additional samples of 2.4 fb^ and 
43.9 fb"^ were collected 30 MeV below the r(35) [below- 
r(35)] and 40 MeV below the r{AS) [below-T(4S')] res- 
onances, respectively and are used for background and 
calibration studies. The trajectories of charged particles 
are reconstructed using a combination of five layers of 
double-sided silicon strip detectors and a 40-layer drift 
chamber, all operated inside the 1.5-T magnetic field of 
a superconducting solenoid. Photons are detected using 
a CsI(Tl) electromagnetic calorimeter (EMC), which is 
also inside the coil. The energy resolution for photons 
varies from 2.9% (at 600 MeV) to 2.5% (at 1400 MeV). 

The signal for T(35) -^ jr]b is extracted from a fit to 
the inclusive photon energy spectrum in the CM frame. 
Any reference to photon energy hereafter will be in the 
CM frame, unless otherwise noted. 

The monochromatic photon from the decay appears 
as a peak on top of a smooth non-peaking background 
from continuum (e^e~ — > qq with q — u,d,s,c) events 
and bottomonium decays. Two other processes pro- 



duce peaks in the photon energy spectrum close to 
the signal region. Double radiative decays T{3S) — > 
7Xw(2P);x6j(2P) ^ 7^1^), J = 0,1,2, produce a 
broad peak centered at 760 MeV due to photons from 
decays of the Xbj{2P) states. The peaks from the three 
Xbj(2P) transitions appear merged due to photon en- 
ergy resolution and the Doppler broadening that arises 
from the motion of the x&j(2P) in the CM frame. This 
Xf,j(2P) photon peak is well separated from the signal 
region of interest (around E^ — 900 MeV). We use the 
peak as a tool to verify the optimization of the selection 
criteria and to determine signal reconstruction efficiencies 
and the absolute photon energy scale. The other process 
leading to a peak near 860 MeV in the photon energy 
spectrum is the radiative production of the T{1S) via 



initial state radiation (ISR) e 



+^- 



liSB.'^i^S). Knowl- 



edge of the magnitude and photon energy line shape of 
this background is crucial in extracting the rjb signal. 

We employ a simple set of selection criteria to suppress 
the backgrounds while retaining a high signal efficiency. 
Decays of the r/b via two gluons, expected to be a large 
component of its decay modes, have high track multi- 
plicity. Hadronic events are selected by requiring four or 
more charged tracks in the event and that the ratio of the 
second to zeroth Fox- Wolfram moments [6[ be less than 
0.98. 

Photon candidates are required to be isolated from all 
charged tracks. To ensure that their shapes are con- 
sistent with an electromagnetic shower, the lateral mo- 
ments [7[ are required to be less than 0.55. The signal 
photon candidate is required to lie in the central angu- 
lar region of the EMC, -0.762 < cos(6'^,lab) < 0.890, 
where Oj^lab is the angle between the photon and the 
beam axis in the laboratory frame. This requirement 
ensures high reconstruction efficiency and good energy 
resolution, and reduces the contributions of ISR photons 
from e+e^ —f jisnXi^S) events. 

Due to the fact that there is no preferred direction in 
the decay of the spin-zero r]b, the correlation of the di- 
rection of the photon momentum in the CM frame with 
the thrust axis [8| of the r]b is small. In contrast, there 
is a strong correlation between the photon direction and 
thrust axis in continuum events. The thrust axis is com- 
puted with all charged tracks and neutral calorimeter 
clusters in the event, with the exception of the signal 



photon candidate. We require |cos0t| < 0.7 to reduce 
continuum background, where 9t is the angle between 
the thrust axis and the signal photon candidate in the 
CM frame. 

Photons from ■n'^ decays are one of the main sources 
of background. A signal photon candidate is rejected if 
it combines with another photon in the event to form a 
TT^ candidate within 15 MeV/c^ of the nominal tt" mass. 
To maintain high signal efficiency, we require the second 
photon of the i:^ candidate to have an energy in the lab- 
oratory frame greater than 50 MeV. 

The above mentioned selection criteria were chosen by 
optimizing the S/ \/{B) ratio between the expected sig- 
nal yield (5) and the background {B). The signal sample 
in the optimization is provided by a detailed Monte Carlo 
(MC) simulation [9]. Since no reliable event generators 
exist to simulate the background photon distribution, es- 
pecially from bottomonium decays, a small fraction (9%) 
of the T{ZS) data is used in the optimization to model 
the background in the region 0.85 < E^ < 0.95 GeV. To 
avoid potential bias, these data are not used in the fi- 
nal fit of the photon energy spectrum. This optimization 
procedure, when applied to the Xbji'^P) yield in data 
in place of the simulated signal, yields the same opti- 
mal selection criteria. The final reconstruction efficiency 
evaluated from the simulated signal MC events is 37%. 

The remaining T{3S) data used for the analysis has an 
integrated luminosity of 25.6 fb~ , which corresponds to 
(109 ± 1) milHon r{3S) events. 

To extract the rjb signal, we perform a binned maxi- 
mum likelihood (ML) fit of the E^ spectrum with 0.5 < 
Ej < 1.1 GeV with four components: non-peaking back- 
ground, Xf)j(2-P) ^ T^Cl-S*), 7/Si?T(lS'), and the r]b sig- 
nal. 

The non-peaking background is parametrized by the 
following probability density function (PDF), f{E-y) = 
A(C-hexp[-a£;^-/3£;2]). 

The form of the Xbj{2P) PDF is complicated by the 
presence of Doppler broadening. Crystal Ball (CB) func- 
tions [10!| are used as phenomenological PDFs for the 
three Xbj{2P) -^ jT{lS) shapes. The CB function is a 
Gaussian modified to have an extended, power-law tail 
on the low (left) side. The relative rates and peak posi- 
tions of the Xf,j(2P) components are fixed to their world- 
averaged (PDG) values [l]|. The parameters describing 
the low-side tail of the CB function are common to all 
three of the XbjC^P) peaks. The Xw(2-P) PDF param- 
eters are determined by fitting the photon energy spec- 
trum, with the signal region (840 to 960 MeV) excluded, 
after subtraction of the non-peaking background. All of 
the Xbj{2P) PDF parameters from this fit, with the ex- 
ception of the overall normalization, are fixed in the ul- 
timate fit to the full photon energy spectrum. 

The PDF of the peaking background from ISR r{lS) 
production is parametrized as a CB function form whose 
parameters are determined from simulated events. To 
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FIG. 1: Inclusive photon energy spectrum in the below-T(4S') 
data, with the non-peaking background subtracted. The peak 
at 1.03 GeV is from the ISR process e+e~ -^ 7/SflT(lS'). 
The superimposed histogram corresponds to a fit with a CB 
function. 



estimate the rate of this continuum component in T(3S) 
data, we use the below-T(3S') and bclow-T(4S') data. 
Figure 1 shows the E^ distribution in the bclow-T(4S') 
data, after applying the selection criteria and subtracting 
the non-peaking background. The fit with a CB func- 
tion yields 35800 ± 1600 events. Extrapolating the cross 
section to the T{3S) energy and correcting for the lu- 
minosity ratio and the small difference in detection ef- 
ficiency, the ISR photon background contribution in the 
final analysis is estimated to be 25200 ± 1700 events. The 
error includes systematic uncertainties. This is consistent 
with and more precise than the rate estimated using the 
below-T(35') data. 

In the final fit of the whole E^ distribution to extract 
the rih signal, all parameters of the Xhj(2P) peak and 
the ISR T{IS) PDFs are fixed to the values from the fits 
described above. 

The rji, signal PDF is a non-rclativistic Brcit-Wigncr 
convolved with a CB function to account for the exper- 
imental Ej resolution. The CB parameters are deter- 
mined from signal MC with the rji, width set to zero. 
Since the width of the rjh is not known, we have chosen a 
nominal value of 10 MeV/c^ for the width. Theoretical 
predictions based on the expected ratio of the two-photon 
and two-gluon widths range from 4 to 20 MeV [l^l ■ The 
free parameters in the fit arc the rji, peak position and sig- 
nal yield, the Xfc./(2P) yield, and all of the non-peaking 
background PDF parameters. 

Figure 2(a) shows the photon energy spectrum and the 
fit result. The non-peaking background is dominant with 
only the prominent Xfcj(2f ) peak visible. In Figure 2(b) 
we show the detail of the signal region, after subtracting 
the non-peaking background. The line shapes of the three 
peaking components, Xbj(2P), ISR T{1S), and the ry^ 
signal are clearly visible. The x^ P^r degree of freedom 
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FIG. 2: (a) Inclusive photon spectrum in the region 0.50 < 
E^ < 1.1 GeV. The component PDFs determined from the fit 
are overlaid on the data points. A prominent Xhj{'2P) peak is 
clearly seen. The dashed line corresponds to the non-peaking 
background component, (b) Inclusive photon spectrum af- 
ter subtracting the non-peaking background, with PDFs for 
X()j(2P) peak (solid), ISR T{1S) (dot), rji, signal (dash) and 
the sum of all three (solid), (c) Inclusive photon spectrum 
after subtracting all components except the rjb signal. The 
CB function shape describes the data points well. 



from the fit is 147/113 = 1.3. Finally Figure 2(c) shows 
the data points with all components except the rji, signal 
subtracted, overlaid with the r]}, signal PDF. The fitted ry^ 
signal yield is 19200 ± 2000 ± 2100 events, where the first 
error is statistical and the second systematic. A total 
systematic uncertainty of 11% is estimated by varying 
the Brcit-Wigncr width in the -qi, PDF to 5, 15, and 20 
MeV, setting the ISR T{1S) component to ±1 a of the 
nominal rate, and varying the PDF parameters fixed in 
the fit by ±1 a. The largest contribution (10%) is from 
the rji, width variation. 

The rji) signal significance is estimated using the ratio 
log(Lniax/-^o)j where imax and Lq are the likelihood val- 
ues obtained from the nominal fit and from a fit with the 
?7f, PDF removed, respectively. Fits have been performed 
where the parameters entering the systematic uncertain- 
ties have been varied within their errors. Data have then 
been fitted with all parameters simultaneously moved by 
one standard deviation in the direction of lower signifi- 
cance. This conservative approach yields a signal signif- 
icance greater than 10 standard deviations. 

As a cross check, we also perform a fit where the yield 
of the ISR T{1S) component is left free, and we obtain 
24800±2300 events for this component. This is consistent 
with the estimate using the below-T(4S') data and pro- 
vides an important validation of the X(,j(2P) line shape 
parameterization. The yield and peak position of the r/b 
signal from this fit are unchanged. 

The E^ signal peak value from the fit is 'dnAt_^\^ MeV. 
We apply a photon energy calibration shift of 3.8 ± 2.0 
MeV, obtained by comparing the fitted position of the 
Xbj{2P) peak to the known PDG value. After including 
an additional systematic uncertainty of 1.3 MeV from 
the fit variations described above, we obtain a value of 
E-^ = 921. 2t^;^ ± 2.4 MeV for the 7jb signal. 

The rjb mass derived from the E^ signal is M{rn,) = 



9388.9I^;3 ± 2.7 MeV/c^ Using the PDG value of 
9460.3 ± 0.3 MeV/c^ for the T(IS') mass, we determine 
the T{lS)-r]b mass splitting to be 71.4l|5±2.7 MeV/c^. 

The value we measure for the splitting is larger than 
most predictions based on potential models [2], but rea- 
sonably in agreement with predictions from lattice calcu- 
lations [13|. The mass splitting between the T{1S) and 
the 77b (IS") is a key ingredient in many theoretical cal- 
culations. The precision of our measurement will allow, 
among others, a more precise determination of the lattice 
spacing [13[ and new precision determinations of as [ij] . 

We estimate the branching fraction by correcting the 
signal yield with the reconstruction efficiency (e) from 
simulated signal MC events, and then dividing it by the 
number of T{3S) events in the data sample. The branch- 
ing fraction of the decay T{3S) -^ jtji, is found to be 
(4.8±0.5± 1.2) X 10"''^, where the first uncertainty is sta- 
tistical and the second systematic. The systematic uncer- 
tainty of 25% comes from uncertainties in the signal yield 
(11%) and e (22%). The latter is obtained by comparing 



the yield of Xbj{'^P) in data to the number of expected 
events, which is calculated from the known branching 



fractions [ll|, the number of T{3S) events, and MC re- 
construction efficiency of X{,j(2P). They show a 13% 
discrepancy, but are consistent within the errors. We as- 
sign the full difference to the systematic uncertainty. A 
total uncertainty in e is obtained, after adding the uncer- 
tainties in the Xftj(2P) branching fractions (18%). 

In conclusion, we have observed, with a significance of 
10 standard deviations, the radiative decay of the T(3iS') 
to a narrow state lying slightly below the T{1S). The 
most likely interpretation of the signal peak is the T{3S) 
transition to the bottomonium ground state, although 
other hypotheses, such as a radiative transition to a light 
Higgs boson, are not excluded. Under the bottomonium 
interpretation, this is the first evidence for the rji, bot- 
tomonium state, the pseudoscalar partner of the T{1S). 
The mass of the r]b is 9388.9^2:3 ± 2.7 MeV/c^ which 
corresponds to a mass splitting between the T{IS) and 
the T]t of 7I.4I3J ± 2.7 MeV/c^ The estimated branch- 
ing fraction of the decay T{3S) —> 7 77b is found to be 
(4.8 ±0.5 ±1.2) X 10-4. 
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